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I. INTRODUCTION 

In special relativity, Lorentz invariance is extended to actual (accelerated) observers in 
a pointwise manner [1]. This locality assumption originates from Newtonian mechanics [2- 
6]. However, as pointed out by Bohr and Rosenfeld [7, 8], field measurements cannot be 
performed instantaneously; indeed, such determinations involve an average over the past 
history of the accelerated observer. The conflict is resolved in nonlocal special relativity, 
where the acceleration of the observer's world line explicitly appears in a universal kernel 
that acts as the weight function for the observer's memory of its past acceleration [9]. 

The principle of equivalence implies a deep connection between inertia and gravitation. 
It is therefore natural to develop a nonlocal generalization of Einstein's general relativity 
(GR). This has been achieved via the teleparallel equivalent of general relativity (GR||), 
which is a gauge theory of the Abelian group of spacetime translations — see [10-15] and 
the references cited therein. In fact, there is some formal resemblance between GR|| and 
electrodynamics. In electrodynamics, the constitutive relations between (E, B) and (D, H) 
could be nonlocal [16, 17]; in a similar way, one can introduce nonlocality into GR|| via 
constitutive kernels. The simplest possibility would involve a scalar constitutive kernel that 
acts as a weight function for a certain average over past events. These are connected to the 
present through the world function [18] in nonlocal general relativity [19-23]. 

Nonlocal general relativity, or equivalently nonlocal GR , is a tetrad theory of gravity, 
where the gravitational potentials are defined via the tetrad field e tl a {x) and x here represents 
an event in spacetime with coordinates x M = (c£ z , x). The tetrad frame is orthonormal, 
namely, 

where g^ u is the metric tensor of the Weitzenbock spacetime and the Minkowski metric 
tensor rj & a is given by diag(— 1, 1, 1, 1) in our convention. In this paper, Greek indices run 
from to 3, while Latin indices run from 1 to 3. The hatted Greek indices a, (3, etc., refer 
to anholonomic tetrad indices, while /i, u, etc., refer to holonomic spacetime indices. We 
use units such that c = 1, unless otherwise specified; moreover, we set k = 8-7rG/c 4 . 

The geometric framework of nonlocal gravity is based on Weitzenbock spacetime in which 
the tetrad frame field is globally teleparallel; that is, the spacetime is a parallelizable mani- 
fold. We can therefore choose a global Cartesian tetrad frame for which, in addition to the 
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curvature of the Weitzenbock spacetime, the corresponding connection vanishes as well. We 
will work with such a tetrad frame and the corresponding gravitational potentials throughout 
this paper. 

The gravitational field strength in nonlocal gravity is given by 

C^ u a = d^e v a — d u e^ a , (2) 

which is the torsion of the Weitzenbock spacetime in our convention. It is interesting to 
compare this equation with the standard expression for the electromagnetic field strength in 
terms of the gauge potential. A particular combination of the field strength and the tetrad 
frame field, designated as the modified field strength, plays an important part in GR|| as well 
as its nonlocal generalization and is given by 

^ a := 2 ~ *-MH + 2e^ a C u ^ . (3) 
The field equations of nonlocal general relativity may then be expressed as 

W = 0, (4) 

d v W v & = (T & » + E^) . (5) 

With our conventional choice of Cartesian tetrad frame, Eq. (4), which is the exact analog 
of the source-free part of Maxwell's equations, is automatically satisfied. We are thus left 
with our main field equation (5), in which the excitation W' v a is related to the modified 
field strength ' & via a "constitutive" relation. In GR||, this constitutive relation is a linear 
proportionality; that is, & = (\/—g / k) ( £! 1v &■ In our nonlocal ansatz, however, this simple 
constitutive relation is extended to include a nonlocal part that involves a causal scalar 
kernel /C, as described in detail in Refs. [19-23]. The field strength C^ 01 , the modified field 
strength (t^ 01 and the excitation 1-L^ v a are all antisymmetric in their first two indices. The 
nonlocal gravity theory clearly specifies the various quantities upon which the kernel could 
depend [20]; in addition, we expect that the nonlocal kernel — as the function weighing past 
events for the present — would vanish for events that are sufficiently distant in time and space. 
Moreover, is the matter energy-momentum tensor, while is the energy-momentum 
tensor of the gravitational field, defined by 

V=g e & » := --^ & {cjn v "p) + cjwp . (6) 
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The connection between this tensor in GR|| and the energy- momentum pseudotensor of the 
gravitational field in GR has been clarified in Ref . [24] . 

Some of the physical consequences of nonlocal gravity have been explored in Refs. [19- 
23]. Nonlocality can act like dark matter; in fact, nonlocal gravity involves a galactic length 
scale Ao = 10 kpc in order to provide a satisfactory explanation for the flat rotation curves 
of spiral galaxies. At present, the observational implications of nonlocality for gravitational 
physics in the solar system are essentially negligible, since the size of the solar system is 
very small compared to Ao [21]. 

The main field equation of nonlocal gravity in its linearized form is presented in the fol- 
lowing section. As in general relativity, the general linear approximation to nonlocal gravity 
can be used to study nonlocal Newtonian gravity as well as gravitational radiation. The non- 
local modification of Poisson's equation of Newtonian gravitation theory has been recently 
studied in detail in Ref. [23]. Therefore, the present paper is devoted to the treatment of 
linearized gravitational waves in nonlocal general relativity. As the treatment of linearized 
gravitational waves in GR is well known, in this paper we concentrate instead on an ex- 
amination of the nonlocal deviations of the theory from the standard general relativistic 
analysis. 

Gravitational radiation damping can explain the steady orbital decay rate of the Hulse- 
Taylor binary pulsar as well as similar relativistic binary systems; indeed, this circumstance 
provides indirect evidence for the existence of gravitational waves within the context of GR 
and GR|| [25-27]. Moreover, nonlocal gravity involves a galactic length scale that is very 
much larger than the orbital size of a relativistic binary pulsar. Therefore, we expect that 
the compatibility of the gravitational radiation damping with GR|| [25-27] would still be true 
in the nonlocal generalization of GRii in much the same way as the influence of nonlocality 
on gravitational physics in the solar system is completely negligible at present [21]. Direct 
searches for gravitational waves continue at present — for a recent review see Ref. [28]. 

The plan of this paper is as follows. In section II, we present the general linear approxima- 
tion of nonlocal GR||. In this case, the Newtonian limit was investigated in detail in Ref. [23] 
and the results are briefly described in section III. The treatment of linearized gravitational 
radiation within the framework of nonlocal gravity depends upon the use of appropriate 
nonlocal kernels. These are derived in section IV and applied to the discussion of generation 
and propagation of gravitational waves in sections V-VII. The length scale Aq associated 
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with nonlocal gravity corresponds to gravitational waves of frequency ~ 10 -12 Hz, while 
present observational possibilities involve frequencies > 10~ 8 Hz. Under these conditions, 
we show that the consequences of the nonlocal theory for linearized gravitational radiation 
turn out to be essentially the same as in GR. 



II. NONLOCAL GRAVITY: GENERAL LINEAR APPROXIMATION 

In the linear approximation, nonlocal GR|| can be treated on a background global 
Minkowski spacetime; that is, the corresponding Weitzenbock spacetime may be regarded 
as a slightly perturbed Minkowski spacetime. The deviation of the Cartesian tetrad field 
from the global background inertial axes may therefore be written as 

e l * = fi* + f/>* ll , e» & = 6£-r«- (7) 

We assume that the absolute magnitudes of the nonzero components of xjj^ are so small 
in comparison to unity that the linear weak-field approximation is valid, in which case the 
distinction between holonomic and anholonomic indices disappears as well. Henceforth, we 
deal with the sixteen gravitational potentials ip^. It is convenient to decompose ip^ into 
its symmetric and antisymmetric components, 

:= 2^ (M „), (j)^ := 2^] . (8) 

Then, 

9ixv = + h llu , ip := T]^^ = -h. (9) 
As in GR, it proves useful to introduce the trace-reversed potentials h^, 

= hfu, - -r)nvh , (10) 

so that h = —h and 

1-1 1 

V> = jV + 2^ + • ( n ) 

In terms of the gravitational potentials ift^, the field strength is 

= <9 M Vw - d v tp afl (12) 
and the modified field strength can be expressed as 

^■iiav = —h u [fj, )tr ] — r] v [fj,h a ] Pt P + -<fitj,a,u + Vv[^a]p P ■ (13) 
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We also note that the Einstein tensor is given by 



G 



d £ a 




(14) 



where □ := rj a Pd a dp. 

The gravitational field equations are given in this case by Eq. (5), where d a 'H fl a v = T t 
in the general linear approximation and 



is our nonlocal ansatz. Thus the linearized gravitational field equations reduce, via Eqs. (14) 
and (15), to the following sixteen equations for the sixteen components of ip^ u , namely, 



where T UfJL is not in general a symmetric tensor, but it is conserved, d fJ T utJ ' = 0, due to the 
fact that CC^cri, is in general antisymmetric in its first two indices. 

The nature of the kernel has been discussed in some detail in Refs. [19-23]. In keeping 
with the general linear approximation under consideration here, we assume a linear kernel, 
in which case /C(x, y) is a universal function K of x — y, /C(x, y) = K[x — y). The causality 
requirement is satisfied once K is nonzero only for situations where x M — is a future 
directed timelike or null vector in Minkowski spacetime. That is, event y must be within 
or on the past light cone of event x. This causality requirement implies that K[x — y) is 
proportional to H(x° — y° — |x — y|), where H(s) is the Heaviside (unit step) function such 
that H(s) = 1 for s > and H(s) = for s < 0. 

The gravitational field is expected to vanish at spatial infinity. We assume that the convo- 
lution kernel vanishes for past events that are infinitely distant in space or time. Moreover, 
it is important to note that only the symmetric part of ip^ u involves propagating fields and 
we assume that there are no incoming gravitational waves. With a causal convolution kernel 
in Eq. (16) and using dK/dx u = —dK/dy a , we find 




(15) 




(16) 




(17) 



where 




(18) 
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and kT U(1 + S^ u is a symmetric tensor. The integrand in Eq. (18) is the sum of four terms; 
however, only the temporal (a = 0) integration makes a contribution to S IMll (x) at y° = 
x° — |x — y| and we thus find 

S,v(x) = J K(|x-y|,x-y)C(x°-|x-y|,y) A, (19) 

where we will henceforth refer to |x — y| in the temporal variables as the retardation. 

The integral relationships that we consider are in general of the Fredholm type [29]; 
however, they turn into Volterra integral relations whenever the kernels are causal. As 
shown in Appendix A, K[x — y) has a unique reciprocal causal convolution kernel R(x — y); 
therefore, Eq. (17) can be written as G^ v = U^, where 

U^ix) = nT„^(x) + S^(x) + J R(x - y)[nT UIM (y) + S^(y)} d A y . (20) 

In this equation, we think of T Vii as the matter source and its convolution with R as the 
effective dark matter source. The source-free equation for the amplitude of the gravitational 
perturbation is then given by 

G^(x) = S^(x) + J R(x - y)S^(y) d 4 y . (21) 

It is a significant feature of this linearized theory of nonlocal gravity that causal con- 
volution kernels have reciprocal causal convolution kernels under physically reasonable 
conditions — see Appendix A. It means that nonlocality in the linearized theory can be given 
a consistent physical interpretation in terms of simulated (fake) dark matter source that 
is the convolution of the actual source of the gravitational field with the reciprocal causal 
kernel. 

We emphasize that in this simplest theory of nonlocal gravity, the kernels K and R 
are scalars. Therefore, once these kernels are known in the rest frame of a gravitational 
system, they can be determined in any other inertial system via Lorentz invariance. In 
the background global Minkowski spacetime, imagine a different inertial frame with inertial 
coordinates and fundamental observers (i.e., those at rest in space) that move uniformly 
with respect to the background, so that x' = Ax + 6, where A is a Lorentz matrix and b 
indicates a constant spacetime translation. Then, 



K(x -y) = K'[A(x - y)] , R(x - y) = R![A{x - y)\ . (22) 
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These considerations are necessary in order to deal with the astrophysics of dark matter in 
moving systems, such as, for instance, the Bullet Cluster [32, 33]. 

Let us conclude this general treatment of linear approximation scheme by discussing the 
gauge freedom of the gravitational potentials. Under an infinitesimal coordinate transfor- 
mation x M t— > x'^ = x M — e At (x), we find that ip^ v i— > if;' = ifj^ + e^^ to linear order in e M . 
Moreover, 

= Kv + e H,v + e W ~ ?7M% > <t>'nv = <t>tiv + <W _ e ^,M ( 23 ) 

and we note in passing that 

t = h - 2e a Q . (24) 

As expected from the electro dynamic analogy, the field strength is invariant under 

a gauge transformation; furthermore, the same holds for the modified field strength £ Mcri ,. 
Thus Eqs. (14)-(21) are all gauge invariant. Further discussion of the gauge freedom of 
gravitational potentials is contained in Appendix B. 

We now consider the physical consequences of the linear approximation scheme, namely, 
the correspondence with Newtonian gravity and linearized gravitational waves. 



III. NONLOCAL NEWTONIAN GRAVITY 



The correspondence of nonlocal gravity with the Newtonian theory of gravitation can 
be established in much the same way as in GR. Indeed, we assume that the gravitational 
potentials ip^ form a diagonal matrix that depends only on the Newtonian gravitational 
potential The resulting modified Poisson's equation for $ is independent of the speed of 
light c; therefore, one may think of the correspondence with Newtonian gravity as a formal 
transition in which c — > oo. 

We start with Eq. (16) and impose the transverse gauge condition h 1 *" jV = (see Appendix 
B). Assuming that the only nonzero component of that is relevant here is hoo = —4:Q/c 2 , 
we find that 2^^ = h^ u is given by (— 2$/c 2 ) diag(l, 1, 1, 1). The only significant component 
of Eq. (16) is thus the \i — v = one with c 2 £o«o — 2$^ and Too = pc 2 , where p is the matter 
density. Moreover, we assume that in the Newtonian limit, 

K(x-y) = 5(x -y°)x^-y), (25) 
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due to the absence of any retardation effects for c — > oo. The reciprocal kernel is then of a 
similar form [21, 23] 

R(x-y) = S(x -y°)q(x-y) (26) 

and Eq. (16) reduces in this case to the interesting form 

V 2 $(x) = 4ttG [p(x) + J g(x - y)p(y)d 3 y] . (27) 

This nonlocal modification of Poisson's equation has been examined in detail in Ref. [23]. 
Nonlocality appears in Eq. (27) as an extra ("dark") matter source whose density is the 
convolution of matter density with the reciprocal Newtonian kernel. That is, nonlocality 
simulates dark matter [34-36]. It is remarkable that, within the phenomenological Tohline- 
Kuhn approach to the problem of rotation curves of spiral galaxies [37-39], an equation 
of the form (27) was suggested by Kuhn in order to account for dark matter as modified 
gravity. In fact, Tohline [37] first suggested that the Newtonian gravitational potential for 
a point mass M be replaced by 

GM GM /|x|\ 
$(x) = --- + — In U , (28) 



where Ao is a constant length of order 1 kpc. Tohline's purely phenomenological proposal 
was later extended by Kuhn and his collaborators in order to resolve the problem of dark 
matter in galaxies and clusters of galaxies [38]. In particular, Kuhn suggested a modification 
of Poisson's equation of the form (27) with the kernel 

Q(x-y) = 1 1 (29) 
47tA |x- y| 2 

This Kuhn kernel Q is such that for p(x) = M<5(x), Eq. (28) is a solution of the modified 
Poisson equation (27). Further implications of the Tohline-Kuhn approach are discussed in 
Bekenstein's lucid review article [39] and, in connection with nonlocal gravity, in Refs. [19- 
23]. In particular, we should mention here that the Tohline-Kuhn scheme is in seeming 
conflict with the Tully-Fisher relation [40]. As discussed in Ref. [23], we take the view that a 
purely gravitational treatment does not contain enough physics to deal with electromagnetic 
radiation aspects needed for a fair comparison of the implications of nonlocal gravity with 
the empirical Tully-Fisher law [40]. 

To determine the kernel of nonlocal general relativity from comparison with observational 
data, it is necessary to extend the Kuhn kernel over all space. This issue has been extensively 
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treated in Ref . [23] ; in particular, two examples were completely worked out explicitly. These 
can be expressed via the reciprocal kernel g(x — y) as 

i i+f^LHO e _ (30) 

47rA (a + uy 

_ 1 l + a(a + u) 

12 — — 7 / — : — s — e , (dlj 

47tA tt(a + w) 

where w = |x — y| and Ao = 10 kpc. Typical values of the parameters a and a are a^ 1 = 
10 Ao and a = 10~ 3 Ao, so that in general aAo and a/Ao are positive constants that are 
small compared to unity and < aa <C 1. We note that q% and q2 are smooth positive 
functions that rapidly fall off to zero at infinity; moreover, they are integrable as well as 
square integrable. The explicit forms of the corresponding \ kernels have been numerically 
determined and presented in Ref. [23]. It follows from the numerical results that xi an d 
X2 drop off with distance u extremely fast and are essentially zero beyond around 2.5 A . 
To deal theoretically with isolated astronomical systems, one routinely assumes that they 
possess sharp boundaries and are therefore compactly supported. In a similar way, we expect 
that in practice the dark counterpart of an isolated astronomical system is more extended 
but still can be cut off beyond a certain distance away from the source and hence considered 
isolated as well. 

The novel features of nonlocal gravity emerge on galactic scales. That is, the deviation 
of the nonlocal theory from general relativity is associated with a length scale of order 
Ao- Indeed, nonlocality disappears for Ao — > oo. Newton's theory of gravitation and its 
relativistic extension in classical GR are devoid of any intrinsic length scale [41]; moreover, 
GR is in good agreement with current solar system data. Whether this situation extends to 
galactic scales is an open problem, however. We assume that what appears as dark matter 
in astrophysics is in fact mainly a manifestation of the nonlocal aspect of the gravitational 
interaction. 

Turning now to the problem of gravitational radiation in nonlocal general relativity, the 
first task before us is the determination of the reciprocal kernel R(x — y) and hence the 
kernel of the theory K(x — y). This is treated in the next section. 
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IV. KERNELS 

In the linear weak-field approximation, we have neglected possible nonlinearities in kernel 
K, and have assumed that our convolution kernel K[x — y) is such that it has a reciprocal 
R(x — y). In the Newtonian regime, it is the reciprocal kernel that can be directly compared 
with observational data. Following our general treatment in Ref. [23], the relation between 
K and R may be expressed — via Fredholm integral equations of the second kind [29] — as 

Q(x) + J K(x - y)g(y)d 4 y = F{x) , (32) 

T(x) + J R(x- y)F(y)d A y = Q{x) . (33) 

The Liouville-Neumann method of successive substitutions makes it formally possible to 
relate these equations [29]. However, we found by way of several trials in Ref. [23] that in 
the Newtonian regime, where the kernels are given by Eqs. (25) and (26), the corresponding 
Neumann series would converge only in situations that were outside the physical domain of 
interest and so we had to resort to the Fourier transform method. 

Beyond the Newtonian limit, as explained in detail in Appendix A, the causality require- 
ment renders K and R Volterra kernels and Eqs. (32) and (33) Volterra integral equations 
of the second kind. The imposition of causality makes it possible to show, via the Liouville- 
Neumann method applied to the Volterra algebra that a unique reciprocal kernel exists and 
is causal [30, 31] — see Appendix A. Restricting our treatment to continuous, absolutely in- 
tegrable and square integrable causal kernels, it is in principle also possible to apply the 
Fourier transform method. 

Let 

/(e) = f /(s)e-«* d'x (34) 
be the Fourier transform of /, where £ • x :— rj a ^ a x^ '. Then, 

f( x ) = 7W-r, f /(Oe** A • (35) 



(2vr) 4 , 

As described in detail in Ref. [23], we work in the space of functions for which such operations 
are permissible. It follows from Eqs. (32) and (33) via Fourier transformation that 

* ({) = __««L, = --*(«- , (36) 

l + fl(£) i + ff£) 
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since the kernels are reciprocal. Starting with kernel R, one can find K provided l+i?(£) ^ 0, 
etc. In Ref. [23], working in the Newtonian regime (c — > oo), we showed that for an 
appropriate reciprocal kernel q, which would be a natural extension of the Kuhn kernel Q 
to all space, kernel \ exists with 

xm = -^mr (37) 

provided 1 + q{\£\) 7^ 0. In particular, for the class of kernels of either form q\ given by 
Eq. (30) or q 2 given by Eq. (31), we showed that [23] 

l + g(|£|)>0. (38) 

To proceed further, we need to find the functional form of R and K in order to go beyond 
the Newtonian limit (c = 00). Therefore, in what follows we wish to concentrate on the 
Dirac delta function that appears in the Newtonian kernels (25) and (26). In general, we 
expect that constitutive kernels R and K should decay exponentially for events that are 
distant in space and time. The spatial exponential decay is already evident in Eqs. (30) 
and (31). To see how a corresponding temporal decay could come about, we recall here a 
mathematical result, namely, that for n — 1, 2, 3, the functions 

S n (s) = ne- ns H(s) (39) 

form a Dirac sequence such that as n — > 00, 5 n (s) — >■ 5(s — + ). It follows that for a positive 
constant length 1/A, we have that as c — >• 00, 

H(t x - t v - i|x- y\)A c e - Ac ^-^ 5{t x - t y ) , (40) 

where the singularity in 5(t x — t y ) is at t x —t y = + . Here the role of n is formally played 
by a fixed positive constant times c. The retardation in Eq. (40) is proportional to 1/c and 
tends to zero as c formally approaches 00. Moreover, we note that 5(t x ) = c 5(x°). Based 
on these considerations, we will henceforth assume that 

R(x -y)= H(x° - y° - |x - y\)A e" A(a;0 - y0) g(x - y) , (41) 

where 1/A is a constant length of order Ao; in fact, for the estimates in this paper we set 
A = a = 1/(10 Ao). Clearly in the limit as c — > 00, we recover Eq. (26). Moreover, we 
assume that g(x — y) is of the type of qi or q 2 for which Eq. (38) is satisfied. It follows from 
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these considerations that R(t x ,x) is a positive function on as well as within the future light 
cone and zero otherwise; moreover, it is integrable as well as square integrable, and rapidly 
decreases to zero at infinity. 

Apropos of the functional form of our scalar convolution kernel (41), it is interesting to 
digress here and mention a solution \1/ of the massless scalar wave equation, n& = 0, such 
that \l/(ct,x) = P(r) exp (—act) with r = |x| and constant a > 0. Then, V 2 P = a 2 P and 
rP(r) is a linear combination of exp (— or) and exp (ar). Thus a possible solution for \1/ is 
of the form r _1 exp [— a(ct + r)], which is in some ways reminiscent of Eq. (41). 

The kernel given in Eq. (41) is causal by construction; its reciprocal, namely, K, is 
causal as well. The proof of this assertion in a rather general context using the Liouville- 
Neumann method is essentially due to M. Riesz [30, 31] and is briefly described in Appendix 
A. The corresponding proof using the general Fourier transform method appears to be 
rather complicated. However, to demonstrate the consistency of our approach in a simple 
mathematical setting that employs only dominant terms as c tends to infinity, let us neglect 
retardation effects and write Eq. (41) instead in the form 

R(x - y) ~ H(x° - y°)A e^^^x - y) . (42) 

It then follows from Fourier transforming Eq. (42) that 

m-j^m), (43) 

which satisfies the requirement that 1 + R ^ 0. Therefore, we can find K using Eq. (36), 

k{Q ^> . (44) 

In calculating the Fourier transform of Eq. (44), it is useful to note that one can use contour 
integration and Jordan's lemma to find 

/OO i^0 x 

Here, in the complex £° plane, there is just a simple pole singularity at —iA(l + q), which 
occurs in the lower half-plane due to the fact that A(l + q) > by assumption — cf. Eq. (38). 
Therefore, 

K(x) ~ -A_H(x°) J q(\Z\y^e- A ^° d s £ . (46) 
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Let us recall that in the limit as c — > oo, we have 

H(x°)A(l + q) e - A ^ x ° . (47) 

Moreover, it follows from Eq. (37) that 

Putting Eqs. (46)-(48) together, we find that, in agreement with Eq. (25), K(x) = 5(x°)x(x) 
in the Newtonian limit, as expected. 

The reciprocity between K and R indicates that our starting point, Eq. (41), is far from 
unique; for instance, Eq. (46) has essentially the same properties as Eq. (42). However, 
we will work with Eqs. (41) and (42) in what follows for the sake of simplicity. Further- 
more, neglecting retardation has led to a manageable expression for the kernel K, namely, 
Eq. (46), which will be employed in the following sections. The nature of this simplification 
is examined in Appendix C. 



V. GRAVITATIONAL WAVES: PROPAGATION 



Let us return to the linearized nonlocal field equations and impose the transverse gauge 
condition \ u = 0. As is well known from GR, this condition simplifies the linearized 
Einstein tensor such that —2G^ V = Dh^ u . The linearized field equations can now be written 

as 

□ V = -2CV(x) , (49) 

where is given by Eq. (20). 

Assuming that the source T^ v is isolated, the corresponding dark source is also then 
expected to be in effect isolated due to the rapid spatial decay of the reciprocal kernel. 
Far from the source and its dark counterpart, the gravitational potentials in the wave zone 
satisfy the field equations 

□ V + + 2 j R(x - y)S, lu (y)d 4 y = . (50) 

Turning now to the expression for S^ u given in Eq. (19), we observe that <£oov = due to the 
antisymmetry of the modified field strength in its first two indices; therefore, Sq^ = 0, but 
S^q is in general nonzero. It follows from Eq. (50) that S 0u = implies that Dh 0u = 0. Far 
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in the wave zone, a fixed detector perceives the emitted gravitational radiation potentials 
ho u to be essentially plane gravitational waves. Therefore, following the same analysis as in 
standard GR, given here in Appendix B for the sake of completeness, one can choose the 
remaining gauge degrees of freedom to set h 0u = and h = 0. Thus with a suitable choice 
of gauge, h^ u is purely spatial and traceless with h l \j = 0. 
Furthermore, it follows from Eq. (50) that 



S l0 (x) + J R(x- y)S t0 (y)d*y = , (51) 
which implies, via Fourier transformation, that Sio(x) = 0, since 1 + R(£) ^ 0. Thus 

S i0 (x) = J K(\x-y\,x-y)£ i ° (x -\x-y\,y)d 3 y = 0, (52) 

where 

C = -\^, j - (53) 

In the decomposition of the gravitational potentials tp^ into symmetric and antisymmetric 
parts, only the symmetric part propagates and the antisymmetric part satisfies Eqs. (52)- 
(53); therefore, to simplify matters, we shall set 



(54) 



for the rest of this paper. 

Summing up, in the wave zone the gravitational potentials reduce to hij(x) such that 

Dhj + 2^(x) + 2Jr(x- y)S lj (y)d 4 y = , (55) 

where in expression (19) for Sij(x), 

since now in Eq. (13) the transverse gauge condition holds; moreover, h 0fl = 0, h = and 
<j>»u = 0. 

To solve Eq. (55), we work in the Fourier domain, where this equation can be expressed 

as 

(to 2 - |k| 2 )%(a;, k) + 2(1 + R)^ = . (57) 

Here u is the wave frequency and k is the wave vector. These form the components of 
a point in the 4D Fourier domain characterized by the propagation vector = (oo/c, k). 
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Using the approach adopted in section IV, where retardation effects are neglected, Sy can 
be calculated from Eqs. (19), (42) and (55), and the result is 

2Sij ~ iAtuqdkDhijiuj, k) . (58) 

Moreover, R is given by Eq. (43). From these results, the dispersion relation for gravitational 
waves takes the approximate form 

u 2 - |k| 2 ~ -iAu(l + R)q(\k\) ■ (59) 

The scale associated with nonlocality, Ao = 10 kpc, corresponds to a characteristic fre- 
quency vq = c/A ~ 10 -12 Hz. We expect that the propagation of such low-frequency grav- 
itational waves would be significantly affected by nonlocality. However, current laboratory 
efforts [28] are directed at detecting gravitational waves of dominant frequency > 1 Hz and 
corresponding dominant wavelength A, where A/Ao < 10~ 12 . Moreover, future space-based 
interferometers may be able to detect low-frequency (~ 10 -4 Hz) gravitational waves from 
astrophysical sources [28]. On the other hand, pulsar timing residuals from an ensemble of 
highly stable pulsars can be used to search for a stochastic background of very low-frequency 
(~ several nHz) gravitational waves [28]. Thus current observational possibilities all involve 
gravitational radiation of wavelengths that are short compared with Ao [28] ; indeed, in all 
cases of interest, 

A < 10- 4 . (60) 

For such wavelengths the contribution of nonlocality to the dispersion relation (59) is 
negligibly small, so that in practice u 2 ~ |k| 2 . To see this, we first note that Eq. (43) 
implies 

l^k)l^ (A2+ ^ 2)1/2 lg(lk[)l, (61) 

where only positive square roots are considered in this paper. Moreover, to simplify matters 
still further, we can estimate |g(|k|)| using the Kuhn kernel instead. For the Kuhn kernel Q 
given by Eq. (29), we have that 

Thus with A = a = 1/(10 A ), A/to = (20tt)- 1 A/A ; hence, we find \R\ ~ (80tt)- 1 (A/A ) 2 , 
which is negligible compared to unity. Similarly, the ratio of the nonlocal term in the 
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dispersion relation (59) to uj 2 has the absolute magnitude « Au; -1 |<7|, which is again of order 
(807r) _1 (A/Ao) 2 ; hence, the nonlocal contribution to Eq. (59) can be neglected. 

Unless observational efforts extend to extremely low frequencies of order 10~ 12 Hz, we 
can safely neglect nonlocal terms in Eq. (55), so that nhij ~ 0. Therefore, for current 
observational possibilities, gravitational wave propagation in the nonlocal theory reduces 
essentially to that of GR and we can recover the TT gauge and hence the two independent 
helicity states of gravitational radiation — see Appendix B. 

We now turn to the influence of the dark matter source upon the amplitude of emitted 
gravitational radiation. 



VI. GRAVITATIONAL WAVES: DARK SOURCE 



The general linear approximation of section II is such that d tl G^ v = in Eq. (14) and 
jvn = o in Eq. (16). Based on the considerations of section V, we can now essentially 
ignore S^ v in Eqs. (17) and (20). This implies that in the linearized field equations given by 
Eq. (49), we have U^ u ~ nT^u, where 7^, consists of the contribution of the source and the 
dark source, namely, 

T^:=T^ + T^. (63) 
Here, the dark source is represented by 

Tj?(x) = J R{x- y)T^{y)d A y (64) 

and both energy-momentum tensors are now symmetric and independently conserved in 
the approximation scheme under consideration here. We will treat the source and its dark 
counterpart as isolated astronomical systems. 

We are interested in the special retarded solution of the linearized field equations given 

by 

tV- v / o n _ K f T^(x° - |x-y|,y) 3 

h (x,x)«— / : : dy. (65) 

in J |x-y| 

Far away from the source, we can introduce in Eq. (65) the approximation that |x — y| 
|x| = r; that is, 

7T(x°, x) » / T< V - r, y) d'y , (66) 
Z7rr J 

so that the solution takes the form of a spherical wave approaching the detector. Further- 
more, let n denote the unit vector that represents the direction of propagation from the 
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source to the receiver. Far in the wave zone, the spherical wave front can be locally ap- 
proximated by a plane wave front with wave vector k = con. The wave amplitude in the 
TT gauge — see Appendix B — can be extracted from Eq. (66) by means of the projection 
operator P l j = 5j — n l rij, namely, 

h% T = (P\P j m - ~P ij P em )h im . (67) 

For instance, if the spatial frame is oriented such that n points in the positive x 3 direction, 
the only nonzero components of h^ T are W^ T = —h^ T = (h — h 22 )/2 and h^ T = h^ T = h U . 

Next, we recall that the conserved symmetric energy-momentum tensor T^ u of an isolated 
system satisfies Laue's theorem, namely, 

J T iJ (ct,x) d 3 x = J r >0 (ct, X )x l x j d 3 x. (68) 

It follows from Eqs. (66)-(68) that, among other things, h^ T will depend upon the second 
temporal derivative of the total quadrupole moment of the system. We note here for the 
sake of completeness that the temporal coordinate of the quadrupole moment is in fact the 
"retarded" time ct — r, where r can be treated as a constant for the purposes of the present 
discussion. Moreover, the quadrupole moment could just as well be replaced by the reduced 
(i.e., traceless) quadrupole moment in the expression for h l ^ T . 

We are particularly interested here in the contribution of the dark quadrupole moment 

Q%{x°) := j T°°(x°, x)xV d 3 x (69) 

to the total quadrupole moment of the system. It follows from the definition of dark energy 
density that 

Q%{x°) = j d A y T 00 (y) j R(x° - y°, x - y)xV d 3 x . (70) 

Let us introduce the new variable u = x — y in Eq. (70); then, using expression (42) for the 
reciprocal kernel and recalling that q(u) is assumed to be spherically symmetric and of the 
general form of either qi or q2, we find 

g^ )~M^' + ^V), (7i) 



where 

M = ! T 00 (v) d 3 v 

M 



M=lT™(y)d 3 y, ^ = / q(u) d 3 u (72) 
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and 

q v ■= J g( u )«V d 3 u. (73) 

Here M is the mass-energy of the radiating system, which is conserved at the linear order, 
and Md is a rough estimate for the corresponding "dark" mass-energy. That is, Md would 
be the dark mass-energy for a pointlike source of mass-energy M; in fact, M^/M « 2/ (aAo) 
for either q\ or q2 with aAo = 0.1 — see Ref. [23]. The time-dependent part of the "dark" 
quadrupole moment may be written as 

£ 4 V) = A f H(x° - y )e- A{x °- y0) Q lj (y ) dy° , (74) 

which is a certain average of the quadrupole moment <2 JJ (t) of the system. 

Suppose we are interested in an astronomical system whose quadrupole moment varies 
with time with a dominant frequency fi that can be detected on Earth via a gravitational 
wave detector within a reasonable span of time; therefore, we expect that fi 3> A. Expressing 
such a Fourier component of Q l ^{y°) as a constant amplitude times cos (fiy° + tp), where tp 
is a constant phase, it follows from 

A J H(x° - y )e~ A( - x °-y^ cos (Vly + <p) dy° = A2 + Q2 [A cos (fix + <p) + Q sin (fix + <p)] 

(75) 

that the relative contribution of the time-dependent part of the dark quadrupole moment 
will be reduced at least by a factor of A/Q 1. We therefore conclude that the contribution 
of the dark source to hji T is essentially negligible for all systems that are currently under 
consideration as possible sources of gravitational radiation that could be detectable in the 
near future. 

VII. GRAVITATIONAL RADIATION FLUX 

Finally, we must compute the flux of gravitational radiation energy at the detector. 
The energy-momentum tensor of gravitational waves S^ u can be calculated in the linear 
approximation using Eq. (6); that is, £^ v is the linearized form of the energy- momentum 
tensor of the gravitational field. In this calculation, we assume, as before, that (f)^ = 0; 
moreover, of the various gauge conditions leading to the TT gauge, at this point only h " )V = 
and h = are explicitly imposed here for the sake of simplicity. Hence, we find that 
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and 



k S, 




(77) 



Let us note that is traceless and that in the absence of nonlocality (i.e., K — 0), it 
becomes symmetric as well, and its form is then reminiscent of the electromagnetic energy- 
momentum tensor, as expected from the linearized GR|| theory. This analogy with electro- 
dynamics regarding the field energy and momentum has been discussed in a general context 
in Ref. [43]; moreover, it has been shown there that in GR||, the energy- momentum tensor 
vanishes for a class of exact plane-fronted gravitational waves. These exact solutions are 
therefore physically meaningless. This lack of physical significance of such exact plane waves 
has no bearing on our local plane waves, as we deal with the circumstance that, far from 
the source, the spherical gravitational waves can be approximated locally — i.e., near the 
receiver — by plane waves. 

We determine the energy flux implied by Eq. (77) in two steps. First, we compare the 
local (K = 0) result Sffi to the Landau-Lifshitz energy-momentum tensor t^ v of gravita- 
tional waves in the linear approximation of GR. Then, we give an estimate of the nonlocal 
contribution to £^ w . 

Under the same explicit gauge conditions for the deviation of the metric tensor from the 
Minkowski metric tensor, namely, h = and hF" \ u = 0, the corresponding Landau-Lifshitz 
tensor tp V , which is in general symmetric, but not traceless, is given by (see the Appendix 



It is straightforward to show from Eq. (77) that S^J , which is defined to be £^ u for K = 0, 
is symmetric and traceless, and can be expressed as 



Each of the expressions in Eqs. (78) and (79) consists of the same three parts, but they 



of Ref. [44]) 




(78) 




(79) 
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differ in the overall numerical factors in front of the first and last parts: These are both ~ 
in Eq. (78), but \ in Eq. (79). 

For the calculation of the energy flux, we impose the additional gauge condition that 
hoy = 0. It then follows from Eqs. (78) and(79) that 

4?-*ofc = ^/W^',o, (80) 

which vanishes in the TT gauge due to the transverse nature of the radiation. For instance, 
if the spatial axes are so oriented that plane waves propagate to the receiver along the a; 3 
axis, then k = 3 in Eq. (80) and h^i = in the TT gauge (see Appendix B). We therefore 
conclude that in the absence of nonlocality, the flux of gravitational radiation energy will 
be the same as in standard GR, in agreement with previous results [26]. 

Let us next consider the contribution of the nonlocal terms in Eq. (77). In the 
Fourier domain, the ratio of the nonlocal term to the corresponding local term is given 
by K = —R/(l + R) in accordance with Eq. (36). However, as explained in section V, 
\R\ ~ (807r) _1 (A/Ao) 2 , which in view of Eq. (60) is completely negligible in comparison to 
unity for radiation of dominant frequency > 10 -8 Hz. We therefore conclude that the non- 
local contribution to the energy-momentum tensor of gravitational radiation can be ignored 
for gravitational waves that may be detectable in the foreseeable future. 

VIII. DISCUSSION 

Experimental efforts are under way to detect gravitational waves in the frequency range 
> 10~ 8 Hz. We have shown that the treatment of such linearized gravitational waves — 
namely, their generation, propagation and detection — in nonlocal general relativity essen- 
tially reduces to that of general relativity. This circumstance is due to the fact that gravi- 
tational waves in such a frequency range have wavelengths that are much smaller than the 
characteristic scale associated with nonlocality, namely, Ao = 10 kpc. For the treatment of 
gravitational waves of frequency < c/Xq ~ 10 -12 Hz, however, nonlocality is expected to 
play a significant role. 
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Appendix A: Causal Kernels 

The purpose of this appendix is to present the chain of arguments that lead us to the 
conclusion that the causal convolution kernel of linearized nonlocal gravity has a reciprocal 
causal convolution kernel. 

1. causality 

The convolution kernel of linearized nonlocal gravity, K(x — y), is causal; that is, it is 
zero unless x — y is a future directed timelike or null vector in Minkowski spacetime. This 
means that x° > y° and 

r]^{x a -y a ){x^-y^<Q. (Al) 

Thus 

{t x -t y ) 2 -\\*-y\ 2 > (A2) 

cr 

and with t x > t y , the causality requirement reduces to 

x°-y°> |x-y| (A3) 

or that 

K(x-y) ocH(x°-y°- |x-y|). (A4) 

The question is whether the reciprocal kernel satisfies a similar relation as Eq. (A4). To show 
this using the Fourier transform method is not simple due to the existence of retardation 
effects — see section IV; however, as discussed below, one can employ the Liouville-Neumann 
method of successive substitutions [29] to demonstrate that under physically reasonable 
conditions the reciprocal kernel is causal as well [30, 31], so that reciprocity holds for causal 
convolution kernels. 

2. Volterra algebra 

The convolution property is independent of causality. To emphasize this point, we start 
in this subsection with a general causal kernel K(x,y). A causal kernel function that is 
continuous on causally ordered sets in Minkowski spacetime is called a Volterra kernel. 
Volterra kernels have some interesting properties that we briefly mention here. 
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We define the product of Volterra kernels K and K' to be 

V(x, y) = J K(x, z)K'(z, y)d i z . (A5) 

The integrand here is nonzero only if 

z°-z°>|x-z|, z° -y° > |z-y| . (A6) 

Summing these two conditions leads to 

x° - y° > |x - z\ + |z - y| > |x - y| (A7) 

by the triangle inequality. Thus V is a Volterra kernel as well. The space of Volterra kernels 
forms an algebra over the causally ordered events in Minkowski spacetime. It is important to 
note that the integration in Eq. (A5) takes place over a domain T>(x, y), which is the region in 
spacetime that is bounded by the past light cone of event x and the future light cone of event 
y, as depicted schematically in Figure 1. For a Volterra kernel K(x,y), it proves useful to 
define iterated (Volterra) kernels K n (x,y) for n = 1,2,3, such that Ki(x,y) = —K(x,y) 
and 

K n+1 (x,y)= [ K n (x,z)K 1 (z,y)d 4 z. (A8) 

JV(x,y) 

It is straightforward to check by a simple change of variable in the integral in Eq. (A5) 
that if K(x, y) and K'(x, y) are convolution kernels and thus just a function of x — y, then V 
is a convolution kernel as well. Moreover, the order of the terms in the integrand of Eq. (A5) 
is immaterial in this case, so that V is also the product of K' and K. Indeed, convolution 
Volterra kernels form a commutative subalgebra of the Volterra algebra. 

Next, consider Volterra kernels K and A; we wish to find a Volterra kernel B such that 

B(x,y)+ / K(x,z)B(z,y)d 4 z = A{x,y) . (A9) 

Jv(x,y) 

This generalized Volterra equation of the second kind has a unique solution in accordance 
with the theorem of M. Riesz [30, 31] given by 

A(x,y)+ R{x ) z)A{z ) y)d 4 z = B(x,y) , (A10) 

Jv{x,y) 

where 

oo 

R(x,y) = Y i K n (x,y). (All) 

n=l 
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FIG. 1: Schematic plot indicating the finite shaded domain T>(x,y) in spacetime. It is the region 
common to the light cone that has its vertex at event x and the light cone that has its vertex at 
event y. 

The Neumann series here converges uniformly on bounded domains and the reciprocal kernel 
R is a Voterra kernel. The proof employs generalized Riemann-Liouville kernels and has been 
sketched in Ref. [31]. The original study of M. Riesz [30] was set in a more general context; 
in our brief presentation here, we have followed the work of Faraut and Viano [31]. 

3. L l and L 2 Volterra convolution kernels 

Our aim now is to reduce the generalized Volterra integral Eqs. (A9) and (A10) to Volterra 
integral Eqs. (32) and (33). To this end, we restrict our attention to Volterra convolution 
kernels that are L l and L 2 functions on spacetime. Indeed, these properties need to hold 
only within and on the future light cone as a consequence of the causality requirement. 
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Let f(x) be any continuous L 1 function over spacetime and define 

JF{x) = jA(x- y)f(y)d 4 y , Q(x) = J B(x - y)f{y)d i y . (A12) 

Then, multiplying Eqs. (A9) and (A10) for L 1 and L 2 convolution kernels by f{y), integrating 
over spacetime and employing Young's inequality for convolutions, we recover Eqs. (32) 
and (33) such that F{x) and Q{x) are continuous L 1 functions over spacetime. Furthermore, 
it is a consequence of Minkowski's integral inequality that if / is L l and A is L 2 , then their 
convolution is L 2 . Therefore, F(x) and Q(x) are L 2 functions over spacetime as well. In 
particular, for our L 1 and L 2 convolution reciprocal kernel we have that 

R(x-y) <xH(x°-y°- |x-y|). (A13) 



Appendix B: Gauge Conditions 

The various gauge conditions imposed in this paper are essentially the same as in standard 
general relativity and their discussion is therefore relegated to this appendix. 
We begin with Eq. (23) and note that in general 

F^F^ + ne". (Bl) 

Let us first find the gauge functions e M for which the transverse gauge condition h ,li )V = 
is satisfied. It follows from Eq. (Bl) that 

n(? = -7T tV . (B2) 

Appropriate solutions of this standard inhomogeneous wave equation can be found to ensure 
that the transverse gauge condition is indeed satisfied. 

Next, we start from such trace- reversed potentials that satisfy h " t „ = 0, but then 
note from Eq. (Bl) that such potentials are not unique. A further gauge transformation 
leads to potentials h! ^ that still satisfy h'^ >u = 0, provided the four gauge functions f^(x) 
satisfy the wave equation 

□ /"(ar) = 0. (B3) 

Let us further assume that Dh 0fl = for plane gravitational waves as in section V. 
We wish to show that the remaining four gauge functions f ^ can now be so chosen as to 
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set h'o^ = and h! = as well. To this end, let us choose the spatial inertial coordinate 
system (x 1 , x 2 , x 3 ) such that the direction of propagation of the waves from the source to the 
detector coincides, for instance, with the positive x 3 direction. This simplifies the analysis 
without any loss in generality. Therefore, near the detector, the wave front can be locally 
approximated by a plane such that h ^ = h ^((), where ( := x 3 — x°. Setting h' ^ = in 
Eq. (23), h! — in Eq. (24) and replacing e M in these equations by /^(C), we find 

^ = h 01 , - = h 02 , (B4) 

d(fo ~ U) -r -r d(f + / 3 ) l-r 

= -^03 = ^00 , -Tz = T, h ■ ( B5 ) 

Let us note that hoo + ^03 = in Eq. (B5) is consistent with the transverse gauge condition 
Ti h ' v = 0, which implies that d(hoo + h 03 )/d£ = 0. This relation can be integrated and 
the integration constant set to zero, as the presence of a nonzero constant here would be 
inconsistent with the fact that these potentials originate from the distant source of gravita- 
tional waves. It is thus evident that can be so chosen as to render h' purely spatial and 
traceless as well. Moreover, we note in connection with the treatment in section VI that 
in this procedure hu — h.22 and h\2 remain invariant; that is, h!\\ — h'22 = hu — h-n and 
tin = hi 2 . 

Finally, let us assume that in addition to the gauge conditions already discussed above, tij 
is such that nhij = 0. Again, near the receiver in the wave zone, the spherical gravitational 
waves associated with these potentials locally behave as plane waves and, as before, we can 
assume that tij = hij((). Then, ti\j = implies that dh^/dC, = and hence ^(C) = 0. In 
this way, we recover the TT gauge of GR, where the two independent states of gravitational 
radiation are given by hu = — /i 2 2 and h±2 = h,2i- 

Appendix C: Fourier Transform of the Reciprocal Kernel 

The purpose of this appendix is to compute 

R{w, k) = J R(t, x)e^- ik x dtd 3 x , (CI) 
where in accordance with Eq. (41), 

R(t, x) = H{t - r)Ae- At q{r) , (C2) 
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A > 0, r = |x| and g is either q± or g 2 given explicitly in section III. Integrating over the 
temporal interval t : r — > oo, we find 

£(w,k) = y g(r)e" ik - x J(u;,r)rf 3 x, (C3) 

where 

J(w, r) := e - (A " iaj)r (C4) 

is such that \J(u,r)\ < 1. Neglecting retardation in Eq. (C2), which is the simple approach 
adopted in section IV, amounts to setting J(u,r) equal to unity in Eq. (C3). 
Using spherical polar coordinates and the standard result that 



e -'|k|r cos o sin g d g = 2 ""^'"'' ' , (C5) 



sin( 


k 


r) 




k 


r 



we find 



4tiA f°° 

R(u, k) = — — — / rq(r) sin(lklr) J (a;, rWr . (C6) 

(A - iw)|k| y 

To proceed further, we assume that u — |k| corresponds to frequencies > 10~ 8 Hz and 
A = a = (10 Ao) -1 , so that A/oj < 10~ 6 . To calculate the integral in Eq. (C6), let us note 
that R for u = |k| can be written as 

R= — f 9 tt / rg(r)e- Ar (l-e 2 ^ r )dr. (C7) 



(u; 2 + iuM) Jo 

Let us now use (ft and g 2 , given respectively by Eqs. (30) and (31), for q and rewrite this 
expression in terms of dimensionless quantities 

X := Ar, := ^ > 10 6 , w:= aa . (C8) 

Then, 

«=-WCTW) lim{w) - lW)(a)l (C9) 

where J 1 - ) is formally obtained from for /3 = and < u> -C 1. Moreover, 

I? > = r X{X + 1+ J ] e-^ x dX, if = f° X + 1 + W e-^-^dX , (CIO) 
J (X + w) 2 J X + w 

depending respectively on whether q± or q 2 is used for q in Eq. (C7). Let us recall here that 
w = 10 for the two examples involving q\ and g 2 that were worked out numerically in 
Ref. [23]. 
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It follows from the Riemann-Lebesgue lemma that 1^' — > as (3 — > oo. Indeed, by using 
exp (2i(3X) = (2i(3)~ 1 d(exp (2i(3X))/dX in and repeated integrations by parts, it is 
possible to express as an asymptotic series in powers of 1/(3 as (3 — > oo. Furthermore, 
can be expressed in terms of the exponential integral function — see, for instance, page 
311 of Ref. [42]. We find 

l[ 0) = ~ - (1 + w)e 2w Ei(-2w) , 4° ) = 1 - e 2w Ei(-2w) . (Cll) 

Thus for < w <C 1, we have 

/{ 0) « -\ - [C + In (2 W )] , J<°> » J<°> + 1 , (C12) 

since for < x <C 1, Ei(— x) ~ C + lnx, where C = 0.577... is the Euler-Mascheroni 
constant — see page 927 of Ref. [42]. More precisely, for w = 10 -4 , l[ ^ ~ 7.443 and ~ 
8.442, while for w = 10~ 3 , /{ 0) w 5.156 and 4° } w 6.151. 

We are now in a position to compare the above exact expression for \R\ in the u = |k| 
case with the estimate used in section V, namely, \R\ ~ 5tt/(3 2 . We find that for [3 as in 
Eq. (C8), our estimate for \R\ is smaller than the more exact value calculated here by a 
factor of about 2 or 3 for w around 10 -3 or 10~ 4 , respectively. 
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